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ABSTRACT 

We present numerical simulations of disc-planet interactions where the planet opens 
a gravitationally unstable gap in an otherwise gravitationally stable disc. In our disc 
models, where the outer gap edge can be unstable to global spiral modes, we find 
that as we increase the surface density scale the gap becomes more unstable and the 
planet migrates outwards more rapidly. We show that the positive torque is provided 
by material brought into the planet's coorbital region by the spiral arms. This material 
is expected to execute horseshoe turns upon approaching the planet and hence torque 
it. Our results suggest that standard type 11 migration, applicable to giant planets in 
non-self-gravitating viscous discs, is likely to be significantly modified in massive discs 
when gravitational instabilities associated with the gap occur. 
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1 INTRODUCTION 

One of the common approximations in studies of disc- 
planet interactions is to adopt a non-self-gravitating 
dis c. Thus, despite over three decades since the work 
of iGoldreich Tremaind (|l979l . I198Q| ). only a relatively 
small numb er of stud i es have ado pted self - gravitating 
disc mod els (Nelson & Benz 2003a, b; Pierens & Hure 2003; 
Baruteau Masset 2008 : Zhang et al. 2008 : Avliffe Bat3 
20101 ). Disc gravity was included in these works in order to 
obtain a self- consistent treatment of planetary migration. 
The discs remain gravitationally stable. 

Recently, iBaruteau et al.l (|201ll ) and I Michael et al.l 

(|201lh performed numerical simulations of planetary migra- 
tion in gravitationally unstable discs. These studies were 
motivated by the need to understa nd the fate of ^iant plan- 
ets formed by disc instability (Du risen et al.|[200 7). Conse- 
quently, these authors consider massive discs which are grav- 
itationally unstable without the presence of a planet because 
the surfa ce density is h igh enough and the disc is sufficiently 
cool rsee lToomreHl9641 ). 

However, there are other types of gravitational insta- 
bilities in discs. Of particular relevance to disc-planet in- 
teractions is gravitational instab ility associated with in- 
ternal disc edges and grooves (jSellwood KahnI Il99ll : 
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IPapaloizou Savoniie|[l99lh . a s it is known that g iant plan- 
ets open annular gaps in discs (|Lin k, PapaloizoJll986t V 

iMeschiari &; LaughlinI ([2008') suggested planetary gaps 
may become gravitationally unstable by analysing the stabil- 
ity of a prescri bed disc profile. Insta bility was explicitly con- 
firmed by . Lin Papaloizoul (|2011al l via linear and nonlinear 
calculations for ^ aps self-consistently opened by a planet. 

& Papaloizou called these instabilities edge modes since 
they are associated with gap edges. 

In this study we explore the consequence of edge modes 
on planetary migration using hydrodynamic simulations. We 
consider the specific configuration of a giant planet residing 
in a gap with unstable gap edges. Our specific aim is to 
understand how edge modes can modify the standard pic- 
ture of planetary migration expected for giant planets. We 
therefore focus on a small set of simulations rather than a 
full parameter survey. 

This paper is organised as follows. We describe our disc- 
planet models and numerical methods in 32] In 33]we discuss 
the expected stability properties of planetary gaps formed 
in our discs. We present migration results in 311 We find 
that the planet migrates outwards as the gap edge becomes 
increasingly unstable. We analyse a case in JS] to identify 
the required source of positive torque and explain how this 
can be attributed to edge modes. We conclude in ^Slwith a 
discussion of implications and limitations of our results. 
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2 MODEL SETUP 

We consider a two-dimensional self- gravitating protoplane- 
tary disc of mass Md orbiting a central star of mass M* . We 
adopt polar coordinates r = (r, 0) centred on the star and a 
non-rotating reference frame. The governing hydrodynamic 
equations are given in lLin Papaloizoul (|2011al V Units are 
such that M* = G = 1, where G is the gravitational con- 
stant. 

The disc occupies r G [n,ro] = [1,25] and its surface 
density E is initialised to 
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where H{r) is the disc semi-thickness defined below, and the 
surface density scale So is chosen by specifying the Keplerian 
Toomre Q parameter at the outer boundary 



where 



is the Keplerian orbital frequency, and 

Ciso — 
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is the sound-speed profile for a locally isothermal disc. We 
set H(r) = hr and fix the aspect-ratio h = 0.05. The disc is 
also characterised by a uniform kinematic viscosity u = 10~^ 
in dimensionless units. The initial azimuthal velocity V(f) is 
set from centrifugal balance with stellar gravity, disc gravity 
and pressure. The initial radial velocity is set to Vr = 3iy/2r. 

We introduce a planet of fixed mass Mp = 2 X 10"^M* 
on a circular orbit at r = 10 = rp(t = 0). This value of Mp 
corresponds to 2 Jupiter masses if M* = Mq , which opens a 
deep gap leading to type II migration in a typic al non-self- 
gravitating viscous disc (|Lin Papaloizoul 1 19861 ). provided 
no instabilities develop. In this standard picture, migration 
follows the viscous evolution of the gap. 



used in disc-planet simulations, leads to mass accumula- 
tion near the planet and could lead to spurious torques 
from within the planet's Hill radius. In a self- gravitating 
disc, it also causes the effective planetary mass M^, which 
is Mp plus disc ma terial gravitationally bound to the planet 
(ICrida et al.ll2009l). to increase with the disc surface density 
scale (|Lin Papaloizoul [201 Ibl ). 

Equation [5] increases the temperature near r^, thereby 
reducing the mass accumulation and effects above. Physi- 
cally, disc material can be expected to heat up as it falls 
into the planet potential and provide a pressure buffer lim- 
iting further mass accumulation. The use of equation [5] in 
the present work is a simplified prescription for accounting 
for this. However, the gaps opened in a disc with Cs = Ciso 
and with equation [5] are similar. The gap widths are iden- 
tical and the gap depth is < 5% deeper in the case with 



2.2 Numerical methods 

The h ydrodynamic e qu ations are integrated using the FARGO 
code (|Massetll2000al lbl: iBaruteau k Massetll2008l l. The disc 
is divided into Nr x iV^ = 1024 x 2048 zones in radius and 
azimuth, spaced logarithmically and uniformly respectively. 
The resolution is approximately 16 cells per H (or 28 cells 
per rh)- The cells are nearly square (Ar/rA(/) = 1.02). An 
open boundary condition is applied at ri a nd a non- refle cting 
boundary condition at ro (as used by Zhan g et al.l ([2008), see 
also GodonI (| 19961 )1. Self- gravity is solved via a 2D Poisson 
integral, in which a softening length Cg = 0.3H is set to 
prevent divergence and approximately accou nt for the disc's 
vertical thickness (Baruteau Massetll2008l V 

The planet is introduced with zero mass at time t — 
20Po, where Po = 27v /Qk{rp{t = 0)) is the Keplerian pe- 
riod at the planet's initial orbital radius. Its mass is then 
increased smoothly over lOPo to its final value Mp. The 
planet is then allowed to respond to disc gravitational forces 
for t > 30Po. A standard fifth order Runge-Kutta scheme 
is used to integrate its equation of motion. The softening 
length for the planet potential is Cp = 0.6H. 



2.1 Equation of state 

The equation of state (EOS) is locally isothermal, so the 
vertically integrated pressure is p = CgE. Before introducing 
the planet, we set Cs = Ciso- When a planet is present, the 
sound-speed Cs is prescribed as 



hrhpdp 



[{hr)y^ + (/lpdp)7/2]2/7 



(5) 



where Q^p 



GMp/dp, dp = y^\r — rp\'^ + is the softened 
distance to the planet, Vp being its vector position and tp is 
the softening length. The parameter hp controls the increase 
in sound-speed near the planet, relative to the Ciso profile 
above, and is fixed to hp = 0.5. The increase in sound-speed 
is Cs/ciso = 1.54, 1.18 at ru and 2r^ away from the planet 
respectively, where ru — (Mp/SM^Y^^rp is the Hill radius. 
Far away from the planet, the sound-speed becomes close to 

Ciso • 

This EOS was proposed by IPeplihski et^ (|2008al l 
in a series of numerical simulations of type HI migration 
in non- self- gravitating discs. The Ciso profile, commonly 



3 GAP STABILITY 

We first consider gap profiles in disc models with Qo = 
1.5, 1.7, 2.0 to assess their stability. These correspond to 
Keplerian Toomre values at the planet's initial radius of 
Qp = 2.77, 3.14, 3.70 and total disc masses of Ma/M^ = 
0.080, 0.071,0.060, respectively. The discs are gravitation- 
ally stable to axisymmetric perturbations. For smooth ra- 
dial profiles, they are also expected to be gravitationally 
stable to non-axisymmetric perturbations near Vp because 
Qp > 1.5. 

The fundamental quantity for stability discussion 
in a structured barotropic disc is the vortensitjQ (e.g. 



^ Strictly speaking, our discs are not barotropic. However, the 
instabilities of interest are associated with an internal structure 
with characteristic thickness if <C r, but the sound-speed varies 
on a global scale and can be approximated as constant, i.e. strictly 
isothermal and hence barotropic. 
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where = r~^d{r^Q^) /dr is square of the epicycle fre- 
quency and Q — v^/r is the angular velocity. For gaps 
opened by a giant planet, the vortensity profile closely fol- 
lows the Toomre Q profile when the latter is calculated with 
Hi instead of Ofc. Specifically, local extrema in 77 and Q occu r 
app roximately at the same rad ii (|Lin k, Papaloizoull2011al V 

iLin Papaloizoul (|2011al l have shown that planetary 
gaps in massive discs are unstable to global edge modes as- 
sociated with max(?7) or max(Q). For the adopted disc mod- 
els, the background Kepler ian Toomre parameter decreases 
with r, i.e. the disc is more self-gravitating at larger radii. 
Thus edge modes are mostly a ssociated with th e outer disc 
(r > Tp), as found bv lLin Pa paloizou (|2011al V 

Fig. [1] compares the relative surface density perturba- 
tion and the Toomre Q parameter for gaps in the above 
disc models. We expect edge modes to have corotation ra- 
dius Tc at a vortensity maximum, or about 2r^ from the 
planet. This corresponds to a local minimum in relative sur- 
face density perturbation just inside the gap. Assuming the 
coorbital region of a giant planet is such that |r — r^l ^ 2.5rh, 
(jArtvmowic j|2003 : IPaardekooper k, Papaloizoul [2009 Vc is 
just within the outer gap edge. This suggests the develop- 
ment of edge modes could bring over-densities into the coor- 
bital region of the planet. 

Edge modes also require coupling to the outer smooth 
disc via self-gravity, which is easier with lower Qo. The pro- 
files here indicate edge modes will develop more readily with 
increasing disc mass and therefore have more significant ef- 
fect on migration as Qo is lowered. 

Differences in gap profiles are also attributable to differ- 
ent effective planet masses, M^. Without carefully choosing 
Mp and the EOS parameter /ip, it is not possible to have 
exactly the same M'^ in discs of varying Sq. M'p typically 
increases with increasing So, which promotes instability be- 
cause gap edges become sharper. However, increasing sur- 
face density generally favours gravitational instability. Thus, 
planetary gaps in discs with decreasing Qo are increasingly 
unstable, though perhaps more so than if M'^ is held fixed. 

Giant planets can also induce fragmentation in self- 
gravitating discs , a phenom enon previo usly reported in SPH 
calculations ( Ar mitage Sz Hanse n 1999.; iLufkin et al.l 1 20041 ) . 
Indeed, for a test run with Qo — 1.3, we found the plane- 
tary wake, as well as the disc, fragments into clumps and 
no annular gap can be identified. In this case migration is 
expected to be strongly affected by clumps, instead of large- 
scale spiral arms which we will focus on below. 



4 MIGRATION AND DISC STRUCTURE 

The above discussion suggests that edge modes could lead to 
non-axisymmetric disturbances inside the gap. Torques may 
then originate from within the planet's coorbital region. Be- 
cause Tc is expected near the outer edge of the coorbital re- 
gion, edge mode over-densities may undergo horseshoe turns 
upon approaching the planet. Furthermore, fiuid elements 
just inside the separatrix will traverse close to the planet 
when executing U-turns, and this may provide significant 
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Figure 1. One dimensional profiles of the azimuthally averaged 
relative surface density perturbation (top) and Toomre Q in the 
outer disc (bottom). Edge modes require a sufficiently peaked 
max(Q). In our disc models, instability is most prominent in r > 
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Figure 2. Migration of a giant planet in massive discs. The Qo = 
1.5 and Qo = 1.7 models develop edge modes throughout, whereas 
the Qo = 2.0 model first develops vortex modes. 

torques. With this in mind, in this section we aim to iden- 
tify the correlation between migration and gap evolution. 

Fig. [2] shows the instantaneous orbital radius of the 
planet in the above disc models. Migration is non-monotonic 
and can be inwards or outwards on short timescales Pq). 
However, with increasing disc mass, outward migration is 
favoured on timescales of a few lO's of orbits. For Qo = 1.5, 
rp increases by 20% in only 70Po • This is distinct from stan- 
dard type II migration, which is inwards and occurs on much 
longer, viscous timescales (tu = r^/v = O(lO^Po) for our 
choice of u). 

Fig. [3] shows the instantaneous disc-on-planet torques 
during the first 20Po after releasing the planet. Instabilities 
develop within this time frame and cause the torque to be 
positive or negative at a given instant. Up to about t = 40Po, 
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Figure 3. Instantaneous disc-on-planet torques. We have made 
this plot more comprehensible by multiplying the torques by a fac- 
tor f = 1 — exp {—\r — rpp/r^). This reduces local contributions 
from the Hill sphere but does not affect the important feature — 
rapid oscillatory torques — or their behaviour as a function of 
Qo- Despite tapering, instabilities still have significant impact on 
disc-planet torques. 

the torques for = 1.5 and Qo = 1.7 both show large and 
rapid oscillations in comparison to the Qo — 2.0 case. We 
shall see that this is related to edge modes developing for 
Qo = 1.5, 1.7 but not for Qo ^ 2.0. 
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Figure 5. Running-time averages of gap properties: the dimen- 
sionless outer gap depth (top) and the gap asymmetry in units of 
Hill radius (bottom). These quantities are defined in N4.2I 



4.1 Disruption of the outer gap edge 

The behaviour of rp{t) correlates with gap structure, partic- 
ularly the outer gap edge. Fig. [H shows the relative sur- 
face density perturbation when instabilities first develop. 
The least massive disc with Qo = 2.0 develops vortices 
rather than global edge modes (cf. Qo = 1.7). This can be 
expecte d since e dge modes require sufficiently strong self- 
gravity |Lin P apaloizou 2011a). 

In weakly or non-self- gravitating structured discs, un- 
stable modes are associa ted with local vortensity minima 
(|Lin k PaDaloizoull2011bl ) and lead to vortex formation. The 
basic state (Fig. [1]) shows that local min(Q), hence min(?7), 
is located exterior to max((3), and is most pronounced for 
Qo = 2.0. We did observe spiral arms to develop in Qo = 2.0 
later on, but these probably resulted from the vortices per- 
turbing the disc, rather than the linear edge mode instabil- 
ity. The important point with Qo — 2.0 is that instability 
leaves the outer gap edge intact and identifiable, unlike for 
more massive discs. 

In Fig. m the Qo = 1.7 case develops a m = 2 — 3 edge 
mode. Protrusion of the spiral arms into the gap makes the 
outer gap edge less well-defined. The surface density in the 
gap is on average higher in > 0^ than in < The edge 
mode spiral inside the gap and just upstream of the planet 
could provide a significant positive torque as the spiral pat- 
tern approaches the planet (from above in the figure). This 
is consistent with the large positive torque around the time 
of the chosen snapshot seen in Fig. O The outer gap edge for 
Qo = 2.0 is not as disrupted and this results in no secular 
increase in for the simulated time for Qo — 2.0. 

In Fig. m large-scale spirals can still be seen for Qo — 
1.5. Weak fragmentation occurs without a collapse into 
clumps. Like Qo — 1.7, the gap is unclean with significant 
disruption to the outer gap edge. It is no longer a clear fea- 
ture as for Qo = 2.0. 



The surface density plots (Fig. |4|, together with the 
torque plot (Fig. [3]), show that edge modes significantly 
modify torques from those expected for standard type II 
migration. Specifically for our models, edge modes are as- 
sociated with the outer gap edge and lead to a secular in- 
crease in orbital radius (Fig. [2]). We shall see below that this 
is because edge modes provide over-densities inside the gap, 
despite the tendency for giant planets to clear the gap of 
material. 

4.2 Gap evolution 

In this section we use the following prescription to examine 
the evolution of the gap structure. We first calculate the 
azimuthally averaged relative surface density perturbation. 

We define the outer gap edge as re, out > Tp such that 
^5](re,out) = 0, and similarly for the inner gap edge re, in < 
rp. The outer gap depth is defined as 5T^ averaged over r G 
[rp, re, out]. We also define the outer gap widths in units of the 
Hill radius, as WoMt — |re,out — r^l/r^, and the inner gap width 
as Win = |''^e,in — rp\/rh- The gap asymmetry is WoMt — w-m. 
Running time-averaged plots of these quantities are shown 
in Fig. [5] 

Edge modes have a gap-filling effect. The case Qo = 1.5, 
for which the planet immediately migrates outwards, has the 
smallest magnitude of gap depth throughout. The Qo = 2.0 
case has the deepest gap and is non-migrating over the simu- 
lation timescale. In the Qo = 1.7 case, the magnitude of gap 
depth decreases relative to that for Qo = 2.0 after about 
t = 80Po. Note that this corresponds to outward migra- 
tion for Qo = 1.7. These observations suggest that material 
brought into the gap by edge modes, is responsible for out- 
ward migration. 
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Figure 4. Logarithmic relative surface density perturbation, log[E/E(t = 0)], showing gap instabihty as a function of Qo- As Qo is 
lowered, the outer gap edge becomes more unstable. This correlates with an increasing tendency for outward migration. 



We have defined gap asymmetry as the distance from 
the planet to the outer gap edge minus the distance between 
the planet and the inner gap edge. It follows that the more 
negative the gap asymmetry is, the closer the planet is to 
the outer gap edge. Fig. [5] shows the planet is located closer 
to the outer gap edge in discs with lower Qo than in discs 
with higher Qo. Furthermore, the non-migrating Qo = 2.0 
case reaches a constant asymmetry, whereas in the outwards- 
migrating cases asymmetry decreases with time. 

Recall that edge mode spiral arms are associated with 
the outer gap edge. The trend above suggests that the planet 
is moving outer gap edge material inwards via horseshoe 
turns (this would also be consistent with shallower gaps 
with decreasing Qo). This provides a positive torque on the 
planet. If on average this effect dominates over sources of 
negative torques, e.g. Lindblad torques or coincidence of 
an edge mode spiral arm with the outer planetary wake 
(|Lin k Papaloizoul lioTTal l. then the planet should on av- 
erage migrate outwards, i.e. closer to the outer gap edge. 



5 A FIDUCIAL CASE 

In order to understand the mechanism by which edge mode 
spiral arms lead to outward migration, we here focus on the 
case Qo = 1.7. This simulation has been extended to t = 
140Po. Fig.[6]and Fig. [71 summarises this simulation in terms 
of the evolution of the planet's orbit and the disc's surface 
density perturbation. The orbit remains fairly circular (e < 
0.06) with no overall change in eccentricity. However, the 
semi-major axis a increases more noticeably than the orbital 
radius rp. 

In Fig. [6l a rapidly increases towards the end of the 
simulation (t > 130Po)- This is because edge modes eventu- 
ally caused the planet to move into the disc lying beyond the 
original outer gap edge. The planet effectively leaves the gap. 



At this point, we found the surface density contrast ahead 
and be hind the planet conforms to outward type III m i- 
gration ("Masset k Papaloizoul 1 20031 : IPeplihski etalll2008bl l . 
This may have resulted from the fact that edge modes in our 
discs supply over-densities ahead of the planet (see below), 
providing the initial condition, or kick, for outward runaway 
migration. 

Once the planet is in type III migration, edge modes 
become irrelevant. This situation differs from migration sus- 
tained by edge mode spirals associated with the gap edge, 
where the planet can still be seen to reside in an annular 
gap. Henceforth we shall focus on the time frame in which 
this configuration holds (t < 130Po)- 

Fig. [71 shows the outer disc (r > rp) is more unsta- 
ble than the inner disc (r < rp). Large-scale, coherent spi- 
rals are maintained throughout the simulation. The gap is 
unclean with material brought into it by the edge modes. 
Over-densities may lie within the planet's coorbital region 
(2.5r/i > r — rp > 0). Such material is expected to execute 
inward horseshoe turns and exert a positive torque on the 
planet. 

The above effect causes the planet to move outwards, 
towards the outer gap edge where the surface density is 
higher. The planet can then interact with that material by 
moving it through inward horseshoe turns, providing further 
positive torque. The torque magnitude should also increase 
with the migration speed. These effects can result in a posi- 
tive feedback akin to that s een in classic type HI migration 
(jMasset k Papaloizoul l2003l ) , which is consistent with the 
faster- than- linear increase in a for t < 130Po. 

This interaction is local and depends primarily on the 
surface density, so the background Toomre Q parameter, 
which decreases with radius on a global scale, is not expected 
to be of direct relevance in this feedback mechanism apart 
from its dependence on the surface density. 
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Figure 6. Orbital evolution of the planet in the Qo = 1.7 disc, 
in terms of Keplerian semi-major axis a (solid) and eccentricity e 
(dotted). These have been calculated assuming Keplerian ellipses 
without accounting for the disc potential. 
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Figure 7. Overall evolution of the Qo = 1.7 case. The logarithmic 
relative surface density perturbation log [E/E(t = 0)] is shown. 

5.1 Surface density asymmetry 

In Fig. [8] we plot two interaction events between an edge 
mode spiral arm and the planet. In the first, we see that the 
spiral wave disturbance extends into within 2rh upstream of 
the planet. This causes a surface density asymmetry ahead 
and behind the planet, the configuration here correspond- 
ing to positive coorbital torque. This is more apparent in 
the second event: over-density builds up just ahead of the 
planet as material undergoes inward horseshoe turns. The 
fluid shocks while executing the U-turn, since giant planets 
induce shocks close to their orbital radii (|Lin Papaloizoul 
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Figure 8. Two examples of how the passage of an edge mode 
spiral arm by the planet can increase its orbital radius. The disc 
model is Qo = 1.7 and the logarithmic relative surface density 
perturbation log [E/E(t = 0)] is shown. 
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Figure 9. Passage of an edge mode spiral arm by the planet in 
the Qo = 1.7 disc model. The logarithmic relative surface den- 
sity perturbation log [E/E(t = 0)] is shown. In this event the or- 
bital radius decreases because the planetary wake is effectively 
enhanced by the edge mode. Nevertheless, like Fig. [S] this plot 
also demonstrates how edge modes can contribute positively to 
the torque by supplying material into the coorbital region for the 
planet to scatter inwards. 



I2OIOI V A stable gap would have been cleared of material by 
a Jovian mass planet and such an over-density ahead of the 
planet would not exist. 

It is important to note that we have chosen the snap- 
shots in Fig. [HI where the orbital radius rp increases, to 
demonstrate how edge modes can provide a positive torque. 
However, not every passage of the spiral arm increases r^, 
as implied by the non-monotonic migration. An example is 
shown in Fig. (9] Since the outer planetary wake is associ- 
ated with negative Lindblad torques, surface density per- 
turbations due to edge modes may enhance it when the two 
overlap. 

The overall outward migration implies positive torques 
produced by edge modes are on average more significant. 
This may not be surprising since the positive torque comes 
from material cross ing r^. This is also the me chanism for 
type III migration (|Masset Papaloizoul l2003l ) which can 
be much faster than migration due to Lindblad torques. 
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Figure 10. Evolution of the amplitude of the m = 2 (solid) and 
m = 3 (dotted) modes in surface density in the outer disc r G 
[^p,^o] in comparison to the total disc-on-planet torque (dashed) 
and the torque contribution from the Hill sphere (dashed-dot). All 
quantities are running-time averages. Fourier modes have been 
normalised by the axisymmetric amplitude and are plotted on a 
logarithmic scale. The upper plot is Qo = 1-7, in which the planet 
migrates outwards. The lower plot is the non-migrating Qo = 2.0 
case. 



5.2 Torques and mode amplitudes 

To emphasise the significance of edge modes on disc-planet 
torques we can compare the fiducial case Qo = 1-7 to 
the non-migrating case Qo = 2.0. Fig. 1101 shows the time- 
averaged evolution of non-axisymmetric modes in the outer 
disc (r > Tp) and torques in these two cases. The amplitude 
of the m*^ Fourier mode is defined as \Cm/Co\, where 



Crr 



/To /'2 



d(\)^(r^ (j)) exp {—imci 



(8) 



and we plot \Cm/Co\ for m = 2, 3. 

The modes have larger amplitudes for Qo — 1.7 and the 
total torque increases with time 0. The figure shows that 
this is attributed to torques from within the Hill radius. By 
contrast, although the Qo — 2.0 case also develops large- 
scale spirals later on, they are of smaller amplitude and do 
not make the Hill torque positive. The total time-averaged 
torque remains fairly constant. 

Fig. [101 along with Fig. [7]— [9l confirms that large-scale 
spirals which extend into the gap can provide significant 
torques. Edge modes naturally fit this description since they 
are associated with vortensity maxima just inside the gap. 



The total torque values are negative, despite overall outward 
migration, because the running time-average is plotted. This av- 
erage includes the phase of planet release, when the disc-planet 
torque is large and negative. 



Because the outer disc is more unstable, the corotation ra- 
dius Tc of edge modes lies beyond Vp so its pattern speed is 
smaller than the planet's rotation. Thus, over-densities asso- 
ciated with corotation approach the planet from upstream, 
but because Vc is actually within the planet's coorbital re- 
gion, the associated fluid elements are expected to execute 
inward horseshoe turns upon approaching the planet. This 
provides a positive torque on the planet, and if the edge 
mode amplitude is large enough, it may reverse the usual 
tendency for inward migration. 

The picture outlined above is consistent with the fact 
that edge modes are have corotation at vortensity maxima. 
A giant planet can induce shocks very close to its orbital 
radius and vortensity is generated as fluid particles execute 
horseshoe turns across the shock. That is, vortensity maxima 
are associated with fluid on horseshoe orbits. So when non- 
axisymmetric disturbances associated with vortensity rings 
develop (i.e. edge modes), the associated over-densities can 
be expected to execute horseshoe turns. 



5.3 Torque distribution 

The analysis above suggests that torques from within the 
gap are responsible for the gradual increase in Vp or a. Fig. 
II II compares torque densities in the fiducial case to the non- 
migrating case. The torques have been averaged over 30Po 
at two time intervals. 

For t G [50, 80] Po radial plots for both cases show a 
positive torque at r^. By t G [80, 110]Po this torque has 
diminished for Qo — 2.0. This is expected for Jovian planets 
as they open a clean gap leaving little material near Vp to 
torque the planet. However, for Qo — 1.7, this torque is 
maintained (and even slightly increased) by the edge modes 
because they bring material into the co-orbital region. Note 
that the planet in the Qo — 1.7 case is migrating outwards 
for t G [80, 110]Po, so migration itself may also contribute 
to sustaining the torque. 

This torque is caused by over-density ahead in compar- 
ison to that behind the planet, i.e. material flowing radially 
inwards across the planet's orbital radius, as seen in Fig. [51 
However, unlike th e single scattering events by vor tices or 
spirals described in Lin & Papaloizoul (|201Q| . l2011al \ which 
dominate most of the migration, here one has to average over 
many spiral-planet interaction events to see the migration. 



6 SUMMARY AND DISCUSSION 

We have performed self- gravitating disc-planet simulations 
to see the effect of large-scale spiral modes associated with 
a gap opened by a giant planet on its migration. In our 
disc models, edge modes lead to outward migration over 
timescales of a few tens of initial orbital periods. This con- 
trasts to s tandard type H inward s migration on viscous 
timescales (|Lin k, Papaloizoul 1 19861 ). This difference demon- 
strates that gap instabilities significantly affect planetary 
migration in massive discs. 

It is important to note that we have specifically con- 
sidered unstable planetary gaps. This di ffers f rom disc mod- 
els em ployed by iBaruteau et al.l (|201lh and iMichael et al.l 
(|201lh , which are at least twice as massive as our fiducial 
case and develop gravitational instabilities without a planet. 
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Figure 11. Torque densities averaged over t G [50,80]Po (solid) 
and over t G [80, 110]Po (dotted). The upper plots are from the 
Qo = 1-7 case, which display outward migration and the lower 
plots are from the non-migrating Qo = 2.0 case. Vertical lines in 
each plot indicate the Hill radius. 

Furthermore, the planet does not open a gap in their models. 
By contrast, gravitational activity in our discs are entirely 
due to the planet opening a gap. The instability then back- 
reacts on planet migration. 

Baruteau et al. (2011) argued that a single giant planet 
in a massive, gravito-turbulent disc effectively undergoes 
type I migration, which resulted in rapid inwards migration 
in their models. This again contrasts to our less massive 
discs, which are not gravito-turbulent and the interaction 
we consider is that between a planet and large-scale spiral 
waves associated with gap edges. 

In the planet's frame, an edge mode spiral, with coro- 
tation at a vortensity maximum exterior to r^, approaches 
the planet from cj) > However, the disturbance extends 
into the planet's coorbital region, so there are associated 
over-densities that execute inward horseshoe turns, which 
provide a positive torque. The resulting outward migration 
means that this positive torque must on average exceed any 
negative torques. 

In order to sustain this positive torque, edge modes 
need to be sustained to supply material to the planet for 
interaction. This in turn requires the existence of vortensity 
maxima, despite the development of edge modes destroying 
them. However, they are easily regenerated by giant planets 
beca use they induce shocks w hich act as a source of vorten- 
sity (|Lin FapaloizoullioToh . 

Planetary migration observed here closely resem- 
bles classic type HI migration because it relies on 
torques from the coorbital region and the accelerated in- 
crease in semi-major axis indicates a positive f eedback 
(jMasset Fapaloizoul [20031 : iFeplinski et al.ll2008bl ). Accel- 
erated migration may be expected because edge modes move 



the planet toward regions of higher surface density (the gap 
edges) so more material enters the planet's co-orbital region. 
However, other outcomes are also conceivable (see ^6.3|) . 

The main differ ence from the type H I migration orig- 
inally described by iMasset Papaloizou is that the flow 
across the planet's orbit is non-smooth, because edge modes 
provide distinct fluid blobs, rather than a continuous flow 
across Vp. Our results above may be interpreted as a discon- 
tinuous runaway type HI migration. 

Also, where type HI migration usually applies to partial 
gaps and therefore Saturn-mass planets, the Jovian mass 
planet used here resides in a much deeper gap. Radial mass 
flux across the gap is still possible, despite the gap-opening 
effect of Jovian planets, because the unstable edge modes 
protrude the gap edge. 

Our understanding of the effect of edge modes is based 
on the geometry of the gap structure. That is, the existence 
of a vortensity maximum (equivalently Q-maximum) close to 
the gap edge, and its association with global spiral modes. 
Thus the mechanism we have identified, that edge modes 
bring material to the planet for interaction, should be a ro- 
bust phenomenon for planetary gaps in massive discs. 

6.1 Comparison to previous work 

Migration here d iffers to scattering by an edge mode spiral 
arm described by Lin Papaloizo in which a sin- 

gle interaction increases the orbital radius of a Saturn mass 
planet by 40% over a timescale of only 4Po , moving it out of 
its gap (see their Fig. 20). By contrast, the 2 Jupiter mass 
planet simulated here migrates outwards by 20% over a few 
lO's of Po, during most of which it remains in a gap, and 
there are multiple encounters with edge mode spiral arms. 

In both cases the torque comes from material crossing 
the planet's orbital radius, but interaction occurs more read- 
ily with the less massive Saturn mass planet since its gap- 
opening effect is weaker. Furthermore, the equation of state 
employ ed in the preliminary simulation of Lin Fapaloizoul 
(|2011al l allows high surface densities near the planet, poten- 
tially increasing the coorbital torque magnitudes. 

Note that if the planet is first allowed to move after edge 
modes develop, then the initial migration direction would 
depend on the relative position of the edg e mode spiral with 
respect to the planet. This was the case in lLin &; Fapaloizoul 
('2011a'), the planet being released when a spiral arm is just 
upstream. 

However, for the models discussed in this paper, edge 
modes do not develop before planet release (t = 30Po). 
This can be seen by the well-defined max((5) , equivalently 
vortensity maximum, in the gap profiles in Fig. [T] The on- 
set of edge mode s would have destroyed this local maximum 
()Lin PaDaloizoull2011al . Fig. 14). 

6.2 Implications 

A consequence of the above is that formation of a clean gap 
is expected to become increasingly difficult in massive discs. 
Planetary gaps correspond to existence of vortensity max- 
ima, but these become unstable to edge modes in massive 
discs and instability tends to fill the gap. Therefore the stan- 
dard picture and formulae for type II migration should not 
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be applied to gap opening planets in massive protoplanetary 
discs. 

We remark that migration is not only relevant to plan- 
ets in protoplanetary discs. Analogous interactions have 
been discu ssed in the contex t of stars in black hole accre- 
tion discs (jKocsis et al.ll201ll : iMcKernan et al]l201lh . Self- 
gravity could be important in outer regions of these discs, 
so our results may also be relevant to these situations. 



6.3 Outstanding issues 

Our focus in this paper was to identify the mechanism by 
which edge modes may torque up the planet. We have thus 
considered one particular disc model and only varied one pa- 
rameter, the surface density scale. Although the simulation 
timescales were sufficient for our purpose, these are short 
compared to disc lifetimes. 

We expect inter-dependencies between planet, disc 
properties and migration for the specific situation of a planet 
interacting with a gravitationally unstable gap, which was 
induced by it in the first place. In our fiducial simulation 
the planet eventually goes into classic type III migration. 
However, if the planet moves out and leaves its gap with a 
slow enough migration speed, it may open a new gap which 
develops edge mode^ and induces further outward migra- 
tion, or even fragmentation (33}. It is also conceivable for 
some disc parameters, the positive torque from edge modes 
are sufficiently small, so that it may be balanced by inwards 
type II migration on longer timescales. Possible long term 
outcomes need to be explored in a parameter study. 

It should be noted that our adopted disc model is biased 
towards outward migration because the Toomre Q decreases 
outwards (approximately like r~^/^) so edge modes associ- 
ated with the outer disc develop preferentially. This is con- 
sistent with the expectation that typical discs become more 
self- gravitating with increasing radius. However, if the disc 
model was such that the inner disc is more unstable, then 
edge modes should promote inward migration. And if inde- 
pendent edge modes of comparable amplitude develop on 
either side of the gap there could be little migration overall. 

Significant torques originate from material close to the 
planet, so numerical treatment of the Hill sphere is a po- 
tential issue. This concern also applies to standard type HI 
migration, but we remark that the adopted equation of state 
was orig inally designed for num erical studies of type HI mi- 
gration (|PeDlihski et al.ll2008al l. Furthermore, the fully self- 
gravitating discs considered here are wha t is required for ac - 
curate simulations of type HI migration (jCrida et al.ll2009l V 
As a check we have also performed lower resolution simula- 
tions which also resulted in outward migration. 

This EOS mimics heating near the planet but is not 
a quantitative model for the true thermodynamics. If this 
EOS under-estimates the heating, then the positive torque 
could be over-estimated and vice versa. On the other hand, 
edge mode corotation resides outside the Hill radius but still 
inside the gap. Thus the numerical treatment within the Hill 



^ Edge modes also require the existence of a wider disc beyond 
the gap edge, so outward migration should become ineffective 
once the planet moves close the disc outer boundary. 



radius does not affect existence of edge modes and their asso- 
ciated over-densities should affect torques originating from 
the coorbital region in the way described in this paper. 
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